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Abstract: We study the process whereby quantum cosmological perturbations be- 
come classical within inflationary cosmology. By setting up a master-equation for- 
mulation we show how quantum coherence for super- Hubble modes can be destroyed 
by its coupling to the environment provided by sub-Hubble modes. We identify what 
features the sub-Hubble environment must have in order to decohere the longer wave- 
lengths, and identify how the onset of decoherence (and how long it takes) depends on 
the properties of the sub-Hubble physics which forms the environment. Our results 
show that the decoherence process is largely insensitive to the details of the coupling 
between the sub- and super-Hubble scales. They also show how locality implies , 
quite generally, that the decohered density matrix at late times is diagonal in the 
field representation (as is implicitly assumed by extant calculations of inflationary 
density perturbations). Our calculations also imply that decoherence can arise even 
for couplings which are as weak as gravitational in strength. 
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1. Introduction and Summary 

Recent precision observations H, |^ of temperature fluctuations within the Cosmic 
Microwave Background (CMB) agree well IQ with the predictions made for them 
by inflationary models, according to which the Universe once underwent an accel- 
erated expansion during its remote past 0. According to the inflationary picture, 
the observed CMB temperature variations are seeded by tiny primordial density 
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perturbations in the much earher universe, perturbations which began as quantum 
fluctuations during the earher inflationary epoch. 

More precisely, within the inflationary picture the primordial density fluctuations 
derive from the quantum fluctuations, of the quantum inflaton fleld, 

0, whose dynamics drives the signature accelerated expansion. The success of the 
inflationary picture requires these quantum fluctuations to be converted into classical 
spatial variations in the energy density whose amplitude varies stochastically as one 
moves from one Hubble patch to another within the Universe during the much later 
recombination epoch. 

At present, the classicalization of primordial perturbations is understood in terms 
of the properties of the quantum state into which the fluctuation flelds evolve (see, for 
example [§,0). It is known that both scalar and tensor density fluctuations evolve 
during inflation into a particular kind of quantum state - one very similar to the 
squeezed state of quantum optics [|, ^. In particular, quantum expectation values 
of products of fields in a highly squeezed state turn out to be identical to stochastic 
averages calculated from a stochastic distribution of classical field configurations, up 
to corrections which vanish in the limit of infinite squeezing 0, |10|, |ll|. Quantum 



autocorrelation functions like {(j){x)(j){y)) become indistinguishable from autocorrela- 
tions for c-number fields ^{x) computed within a classical stochastic process, when 
computed in the limit of large squeezing. 

Our purpose in this note is to describe what these analyses do not address: the 
decoherence process through which the initial quantum state evolves into the corre- 
sponding stochastic ensemble of field configurations. That is, we wish to understand 
how the density matrix, p[(p{x), ^'{y)] = {ip{x)\p\ip'{y)) , for each of the various quan- 
tum fields, evolves from an initially pure state, 

p[cpix),cf'iy)] = ncp{xW[^'iy)], (1.1) 

describing the initial conditions (where \E'[(/9(x)] = {ip{x)\ip)), to a mixed state, 

p[y.(x), ^'(y)] = P[^ix)] 5[cpix) - cp'iy)] , (1.2) 

describing the later classical stochastic probability distribution for the field ampli- 
tudes. Such a transition is implicit in the standard predictions of infiationary impli- 
cations for the CMB, and our goal is to understand to what extent these predictions 
might depend on the details of how this decoherence is achieved. 

Decoherence is a much-studied process outside of cosmology, although not all 
of the conceptual issues have been resolved. In particular, the transition from pure 
to mixed state we seek typically arises whenever the degrees of freedom of interest 
interact with an 'environment' consisting of other degrees of freedom whose properties 
are not measured. The measured degrees of freedom can make a transition from a 
pure to a mixed state once a partial trace is taken over this environmental sector. 



2 



We are most interested in those modes which are responsible for the correlations 
which are observed in the CMB and in structure formation - modes which have only 
come within the Hubble scale relatively recently. We take the relevant decohering en- 
vironment to consist of those modes having much shorter wavelength, which entered 
the Hubble scale much earlier and whose correlations are not directly visible in the 
sky. By borrowing techniques from atomic and condensed-matter physics we set up 
a master equation for the density matrix of the observed modes which allows us to 
robustly identify how these modes decohere provided only that (z) the environment is 
large enough to not be overly perturbed by the modes whose decoherence we follow; 
(ii) the interactions with the environment are weak; and (ii) the characteristic cor- 
relation time of fluctuations in the environment is much smaller than the timescales 
of interest to the decoherence process. 

We use this formalism to identify how the decoherence depends on the proper- 
ties of the short-wavelength modes which are assumed to make up the environment. 
Other authors have also examined the decoherence problem, often using the formal- 



ism of influence functional [0 (see, for example |]T0|, |T3|, |T5|, |T6|, [T^)- Unfortu- 
nately, the complexity of this formulation often forces calculations to be performed 
only for free fields, which couple to one another only by mixing in their mass or 
kinetic terms, leaving the uncertainty as to whether the results are restricted to the 
special choices which are made for calculational purposes. The main virtue of our 
treatment in terms of the master equation is that it allows us to identify the nature of 
the decoherence process in a controllable approximation which can encompass very 
general, realistic environment-system interactions. 
Our analysis leads to the following results: 

• We find that the quantity in the environment which controls the decoherence 
process is the autocorrelation function of the interaction Hamiltonian, V{t), 
which couples the observed modes to the environment: Tt env[pSV(t) 6V(t')], 
where the trace is over the environmental sector, p denotes the system's den- 
sity matrix and 6V(t) = V(t) — Tronv[p^(^)] represents the fluctuations of 
V{t) (in the interaction representation). More complicated dependence on the 
properties of the environment only arise at higher order in V{t). 

• We find that the generic situation is that the observed modes decohere into 
a density matrix of the form of eq. ( |1.2| ), which is diagonal when expressed 
in terms of field eigenstates. This basis typically plays a special role because 
locality usually dictates that the important interactions between the modes of 
interest and their environment are diagonal when written in this basis. This 
point has previously been emphasized in ref. [O] 



We show that the final stochastic probability distribution in eq. ( |1.2| ) is given 
in terms of the initial pure state by P[lp{x)] = \{{p{x)\ilj)\'^, and thereby jus- 
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tify more precisely the standard arguments which imphcitly use the classical 
stochastic field distribution which reproduces the computed initial quantum 
fluctuations of the inflaton. 

• We find that short wavelength modes of the environment have the potential 
to decohere longer-wavelength modes, but need not do so. In particular they 
should not do so if the environment is prepared with its modes in their adia- 
batic vacua. It is this observation which explains why we are able to measure 
quantum states at all in the lab, given that there are always an enormous num- 
ber of short-wavelength modes to which any given experiment does not have 
access. 

• Wc find that for simple choices for the short-distance environment there is 
ample time for decohcrcncc to occur during or after infiation, even if the modes 
of interest are coupled to the environment with gravitational strength. Since 
all modes couple gravitationally, this shows that decoherence is ubiquitous 
for modes produced during inflation which are now observable on the largest 
scales. In the special case that decoherence arises due to Planck-suppressed 
couplings between super-Hubble inflaton modes and thermal sub-Hubble modes 
after reheating, decoherence is only now beginning to occur for those modes 
which cross the horizon during the recombination epoch, raising the intriguing 
possibility that this might have observable implications. 

We next present these results in more detail. We do so by flrst briefly reviewing 
our approach to decoherence with an emphasis on the limits of validity of our cal- 
culations. We then examine two candidate environments, both of which are known 
to exist in most inflationary scenarios: (i) the short-wavelength modes of the metric 
and inflaton themselves, during inflation and afterwards (for which we flnd our cal- 
culations break down, although wc give arguments as to why decoherence from this 
source may be unlikely); and (ii) those short-wavelength modes which are the flrst to 
thermalize after inflation ends (for which we argue our calculations apply, and show 
that decoherence could reasonably occur in the time available) . 

Our purpose is not to argue that either of these sources must provide the dom- 
inant source of decoherence, but rather to argue that our calculational tools suffice 
to identify that environments exist which have sufficient time to decohere initial 
quantum fiuctuations. Our hope is that these tools can be used to compare po- 
tential sources of decoherence, with a view to identifying those which are the most 
important and how the detailed features of decoherence depend on the nature of the 
environment assumed. In particular, we hope to find observational consequences of 
the decoherence process which might be possible to use to test further the inflation- 
ary paradigm for generating primordial density fluctuations. We group some of the 
details of the calculations into appendices. 
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2. Decoherence and Environments 



The word 'decoherence' can mean a variety of things (see, for example, [T^, ^ 
for extensive discussions), but for the present purposes we intend it to refer to the 
transition from a state which is initially pure — i. e. its density matrix can be written 
P = for some state- vector {ip) — to a mixed state (for which no such state- 

vector exists). Our interest is in how long this transition takes, and in the form of the 
final mixed state to which the system is led. Once the final mixed state is obtained, 
its interpretation as a classical statistical ensemble is most simply found using the 
basis in which it is diagonal, since in this basis we have pred = TrenvP = X]nP'il'^)('^l' 
with < p„ < 1 denoting the classical probability of finding the system in state \n) 
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Convenient diagnostics for determining the purity of the state represented by 
a density matrix p are given by the basis-independent expressions: = p and 
Tr [p^] = 1, either of which are satisfied if and only if the state described by p 
is pure. The equivalence of these criteria for purity rely on the fact that density 
matrices are defined to be hermitian, positive semi-definite and satisfy Tr [p] = 1. 

2.1 The Master Equation 

We now set up a master equation which allows us to describe the pure-to-mixed 
transition quantitatively, and in a way which can be applied to a broad variety of 
choices for both the observable sector and the environment. Our presentation here 



follows closely that of ref. ||22| (see also [jT2|, |23| for alternative approaches). To do 



so we write the system's Hilbert space to be the direct product of the observable 
sector, y4, and the environment, B: Ti. = Ti^ ®'Hb, and write the total Hamiltonian 
governing the system as the sum H = Ha + Hb + V , where Ha and Hb describe the 
dynamics of the subspaces A and B in the absence of their mutual interactions, and 
V is the interaction term which couples them. 

For the applications of present interest we may also imagine the system starts off 
without any correlations between these two sectors, and so is described by an initial 
density matrix: p(t = tg) = Qa® Qb- The presence of the interaction V generates 
correlations between the two sectors as the system evolves; these correlations make 
a general description of their further evolution difficult. A great simplification is 



possible, however, if three conditions are satisfied |24, 25, E6l: 



1. The interactions between sectors A and B are weak; 

2. System B is large enough not to be appreciably perturbed by its interactions 
with system A; 

3. The correlation time, r, over which the autocorrelations {5V{t) 6V{t — T))B are 
appreciably nonzero for the fiuctuations of V in sector B, are sufficiently short. 
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In this case the dynamics of system A over times t ^ t does not retain any memory of 
the correlations with B since these only survive for much shorter times. This allows 
us to treat the evolution of A in the presence of i? as a Markov process. What emerges 
is a picture wherein the interactions with B simply provide a stochastic component 
to the evolution of sector A. From this picture, we can derive a master equation 
which describes how the environment B affects the evolution of the density matrix 
describing sector A, and in particular how quickly sector B leads to decoherence in 
system A. 

It is convenient to work within the interaction representation, for which the time 
evolution of operators is governed by Ha+Hb while the evolution of states is governed 
by V. We concentrate on obtaining an evolution equation for the reduced density 
matrix, pA{t) = Tr B[p(t)], since this includes all of the information concerning mea- 
surements involving observables only in sector A. As Appendix ^ shows, writing the 
interaction Hamiltonian in terms of product basis of operators — V = Aj CS) Bj 
— and using the three assumptions given above allows the derivation of the following 
master equation for pA(t): 

{Bj)B - ^ Wjk[ (^pAAjAk + AjAkPA - 2AkPAA^ | . (2.1) 

where there is an implied sum over any repeated indices, and (• ■ ■ )b = Tr £[(• ■ ■ )qb\ 
denotes the average over the initial (and unchanging) configuration of sector B. The 
assumption of a short correlation time, r, has been used here to write 

{5B^{t)5Bk{t'))B = W,fc(t)5(t - t') , (2.2) 

which defines the 0{t) coefficients W*^ = W^j. 

The crucial observation is that the three assumptions listed above ensure that 
corrections to this equation are of order (V^t)^, where r denotes the correlation time 
over which quantities like {5V{t) 5V{t + t))b are appreciably nonzero. The use of 
perturbation theory to evaluate dpA/dt therefore requires only that r be sufficiently 
short: Vr ^ 1. However once eq. ( p.lj) is justified in this way, its solutions can be 
trusted even if they are integrated over times t for which Vt^ 1. 



OpA 

at 



Pa, Aj 



2.2 Local Interactions and Spatially-Small Fluctuations 

We now specialize to the case of interest for field theory, where interactions are local. 
Suppose, then, that sector A and B interact through interactions of the local form: 

V{t) = J d^x Ai{x,t)Bi{x,t) , (2.3) 

where there is an implied sum on i, Ai denotes a local functional of the fields de- 
scribing the A degrees of freedom, and Bi plays a similar role for sector B. 
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In the special case that the degrees of freedom in sector B have very short corre- 
lation lengths, £, (as well as times, r), their influence on dpA/dt can be represented 
in terms of local 'effective' interactions. That is, suppose 

{6Bi{x, t) 6Bj{x', t'))B = Uij{x, t) 6^{x - x') 6{t - t') , (2.4) 

where U*j = lAji is a calculable local function of position which is of order I'^t. As is 
shown in Appendix]^, under this assumption the master equation, eq. ( p.l|) becomes: 



dpA 
dt 



J d^x pA,Aj {Bj)B-^ J d^xUjk AjAkpA+pAAjAk-2AkpAAj . (2.5) 



Notice that this equation trivially implies dpA/dt = for any pA which commutes 
with all of the Aj. Eqs. ( |2.1|) and ( p.5| ) are the main results in this section, whose 
implications we explore in detail throughout the remainder of the paper. 

We remark that eqs. ( p.2[ ) and ( p.4| ) are a fairly strong conditions, inasmuch as 
they exclude significant correlations along the entire light cone. In particular, they 
exclude the application of the master equations, ( |2.1| ) and (^.5]), to situations where 
{dBj{x,t)6Bk{x' ,t')) depends only on the proper separation, s{x,t]x',t'), between 
(x, t) and {x', t'), such as when B is prepared in a Lorentz-invariant vacuum (like the 
Bunch-Davies vacuum of de Sitter space). 

2.3 Solutions 

If all of the operators Ai{x,t) commute at equal times then the implications of 
eq. ( [2. 51 ) are most easily displayed by taking its matrix elements in a basis of states 
\a)t for which the operators Ai{x,t) are diagonal: Ai{x,t)\a)t = ai{x)\a)t (for fixed 
t). With this choice we denote 

Pt[a,a] = t{a\pA\a)t, (2.6) 
and so using the interaction-picture evolution dt\a)t = —iHA\(y)t, we find 

^[a,a] = {dt{a\)pA\a) + {a\dtpA\a) + {a\pAidt\a)) 

= {a\dtPA\a) + i{a\[pA, HA]\a) . (2.7) 

Notice that the last term in this equation combines with the term involving (V) b in 
eq. (|2.5| ) to give the combination i{a\[pA, Ha + {V) b]\c() in the evolution of the matrix 
element p([Q;,a]. Because these terms describe a Hamiltonian evolution they cannot 
generate a mixed state from one which is initially pure. For our present purposes 



of understanding decoherence we need follow only those terms of eq. ( p.7|) which are 
O(y^) and which do not simply describe pure-state evolution. 
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The general solution to this equation is easily given in those circumstances for 
which the term involving Ha is negligible, since in this case evaluating the matrix 
elements gives 



-^ = -^p^jdz 



adz] - aAz] 



{Bi{z,t))B 



-pt / dh 



adz - adz 



aj[z) — aj[z) 



These equations are readily integrated to yield 



pt[a,a\ 



Pto[a,a]e ^ , 



Uij{z,t) 



(2. 



where 



E 

r 



dt'd^z 



to 



dt'd^z 



to 



adz] - adz] 



adz) - adz 



{B^{z,t'))B 



aj[z) — aj[z) 



(2.9) 



Kijiz,t'). 



This solution describes the evolution of the reduced density matrix over timescales 
for which the unperturbed evolution due to Ha is negligible. 

As expected, the mean-field term (involving (Bi) b) only changes the phase of the 
initial pure state and so it may be neglected to the extent that only the transition 
from pure state to mixed is of interest. The same would be true for the Ha terms in 
eq. ( p.7| ) if [if^j^i] = 0. Otherwise Ha + {V)b acts to evolve, but not to decohere, 
the initial pure state. By contrast, those terms involving the fluctuations Uij cause 
the reduced density matrix to take the form of a classical Gaussian distribution in 
the {ai{x)}, whose time-dependent width is controlled by the local autocorrelation 
function Uij. Provided this width shrinks in time, at late times the system evolves to- 
wards a density matrix which is diagonal in the |a) basis, with diagonal probabilities 
that are time-independent and set by the initial wave-function: 



Pt[a] 



Pt[a,a\ = pto[Oi,a\ 



to 



a 



(2.10) 



Here "^to i'^] is the wave-functional for the initial pure state. 

These solutions describe the decoherence of the initial state into the classical 
stochastic ensemble for the variables, {oj}, which diagonalize the interactions with 
the decohering environment. Since general considerations of locality typically re- 
quire the interactions, Ai to be ultra-local polynomials of the fields, this shows that 
the classical distribution to which decoherence leads is generically an ensemble in 
field space, ip{x) [jT8|, just as is implicitly assumed in standard calculations of 
inflationary perturbations. 
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3. Decoherence in Cosmology 



We next apply the above formalism to inflationary cosmology, in order to identify 
potential sources of decoherence for those fluctuations whose influence is observed in 
CMB measurements. 



3.1 Generalization to FRW Geometries 



We start with a short discussion of how eq. ( p.5| ) must be modified when applied to 
an expanding Friedmann, Robertson Walker (FRW) universe. We may do so either 
using co-moving time or conformal time, leading to physically equivalent results. 

The starting point is the generalization of eq. ( |2.3D to curved space. If we take 
the interactions Ai and Bi to be Lorentz scalars, the insertion of a factor of the 
3-metric's volume element gives 



V{t) = a^{t) / £x Ai{x,t)Bi{x,t) . 



(3.1) 



Keeping in mind that the delta function transforms as a tensor density, eq. (|2] 
becomes 



{SB,{x,t) SB,{x',t'))B 



lAiji^x, t) 



(5^(x - x') 6{t - f) 



where U{x,t) is also a Lorentz scalar. With these choices eq. (|2.5|) generalizes to: 
dpA 



(3.2) 



dt 



d^x 
(fx Ui 



PA,Ai. {Bi)B 

AjAkPA + PAAjAk - 2AkPAAj 



jk 



(3.3) 



The equivalent result in conformal time is obtained using the replacement dt = a dt]. 
3.2 Short- Wavelength Modes as Environment 

For applications to inflationary predictions for the CMB we take sector A to be 
comprised of those modes of the inflaton-metric system whose wavelength allows 
them to be visible within current measurements. In particular, we take A to consist 
of the Mukhanov field PDI, |5D|, 



V = a 



+ —'lb 



(3.4) 



which describes the physically relevant combination of the inflaton fluctuation, 50, 
and the scalar metric fluctuation, ip (see Appendix for more specific definitions). 
Here a{r]) denotes the background scale factor as a function of conformal time, t], 
for which 7i = a' /a, and ip is the background inflaton configuration. Primes denote 
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differentiation with respect to rj. What makes this variable convenient is that its 
hnearized evolution is governed by a ffat-space, free-field action in the presence of a 
time-dependent mass, = —z"/z: 



S2 = l I d^x 



v'^ -Y'divdjV+ ( — )v^ 



(3.5) 



where z = aip' /Ti and is the spatial FRW metric (which we choose in what follows 
to be flat). 

We must then make a choice for the decohering environment, sector B. A min- 
imal choice for B consists of the short-wavelength modes of the infiaton and metric 
fields themselves, since these modes and their mutual interactions are in any case re- 
quired to exist by the assumed fluctuation-generating mechanism. In this section we 
explore the extent to which these interactions can suffice to decohere longer-distance 
infiaton fluctuations, given their very weak strength. 

In particular we ask how the time-scale for decoherence depends on the physical 
parameters of this environment, in order to determine if sufficient time is available for 
decoherence to have been completed for those modes relevant to the CMB between 
horizon exit and re-entry. Our goal in so doing is not to argue that these are the 
most important sources of cosmological decoherence, but rather to show that the 
infiaton couples with sufficient strength in the early universe to plausibly have had 
enough time to have decohered before horizon re-entry. 

Causality Issues 

We have seen that the master equation derived above only applies if the environment's 
correlation time, r, is much smaller than the times over which the system of interest 
evolves. Since the Hubble scale, H, sets the natural timescale for evolution in an 
expanding universe, this formalism is most naturally applied to situations where it 
is sub-Hubble modes {k/a 3> H) which decohere super-Hubble modes {k/a < H). 
Can this type of decoherence be consistent with causality? 



At first sight one might think not, since there are general arguments ||3^ which 
restrict how much short-wavelength quantum fluctuations can modify background 
fields like (0) = yj. However our interest in this paper is in how the system's density 
matrix evolves between horizon exit and re-entry, and in particular how the transition 
occurs from a pure state to one that is mixed. But this is to do with the entanglement 
of short- and long-wavelength modes and the destruction of the phase coherence 
amongst the long-wavelength modes which arises once the short-distance modes are 
traced out. Since mean quantities like tp are not changed by this decoherence, it is not 



constrained by considerations such as those of |3^ . We believe our situation is closer 
to Einstein- Podolsky- Rosen experiments, for which entanglements between widely- 
separated electrons can give rise to seemingly nonlocal evolution of the system's 
density matrix. 
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Figure 1: A cartoon of the evolution of the wavelength of a perturbation and the Hubble 
length vs time during inflation and afterwards. The vertical lines indicate the end of 
inflation, and the beginning of the subsequent reheating epoch. 



3.3 Decoherence After Reheating 

We now must decide when to look for decoherence between the epoch when the 
observed modes leave the horizon during inflation, and their re-entry in the not-so- 
distant past (see Figure p. We seek an instance when the short- wavelength modes 
are prepared in states for which the short- distance assumptions, eqs. (|2.2|) and (|2.4| ), 
are satisfied. 

The simplest such instance occurs after the reheating epoch. Reheating occurs 
once inflation ends and modes begin to re-enter the Hubble scale, since they can 
then acquire more complicated dynamics. Detailed studies indicate that under some 
circumstances this dynamics can lead to the reheating [31| and/or preheating |^| of 



the universe, giving rise to the thermal state which late-time Big Bang cosmology 
assumes. 

Since thermal fluctuations can often satisfy the short- correlation assumptions, 
eqs. (|2.2| ) and (|2.4| ), we now examine whether a thermal bath of sub- Hubble modes 



can decohere the longer-wavelength, super- Hubble modes through their weak mutual 
self-interactions. Within this picture the first modes to thermalize would do so as 
in standard treatments, about which we have nothing new to say. Our interest is 
instead in how this nascent thermal (or other) mixed state acts to decohere those 
modes which only enter the Hubble scale much later, whose imprint we now see in 
the CMB. Since this type of decoherence cannot begin until after inflation ends, a key 
question is whether there is sufficient time over which to decohere the super-Hubble 
inflaton modes before they re-enter. 

Decoherence Rates 

In order to estimate whether there is sufficient time to decohere super-Hubble modes 
we require a simple model for the lowest-dimension inflaton/heat-bath interaction, 
which we take to have the form of eq. (p^.3|), with A = gv and where S is a local 
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operator having engineering dimension (mass)'^, with d > 0. On dimensional grounds 
the couphng g then has dimension (mass)^"*^. With these assumptions the function 
I4{x,t) defined by eq. ( p. 21) has dimension (mass)^'^"^. For later numerical estimates 
there are two cases of interest: (i) a dimension-4 (marginal) infiaton-bath coupling, 
for which d = 3 and g is dimensionless; and (ii) a dimension-5 interaction with d = 4 
and coupling g = 1/M, for some large mass scale M. 

Suppose now that the physics of sector B is homogeneous in space and is char- 
acterized by a single mass scale, A(t), which can be slowly evolving with time as 
the universe expands. For instance, A{t) might be given by the temperature T{t) if 
B were described by a simple thermal state. Working in fiat space we see that on 
dimensional grounds we can take 



{B{x, t))B = ci A^it) and U{x, t) = ca A^'^-^(t) 



(3.6) 



where ci and C2 are calculable dimensionless real numbers (which might also depend 
logarithmically on A, for large A). We imagine C2 > 0, which should be true for a 
broad class of choices for the environment (sector B). 
Using these expressions, we then find 



{V{t))B = gci A'^(t) a%t) / d^x v{x, t) 



and 



dpA 
dt 



{V)b,Pa 

-g''c2A^^-\t)a\t) I d^'x v{x), v{x),pa 



(3.7) 



(3. 



Evaluating this last equation within field eigenstates, \V{x))t, and integrating 
over times short compared with allows the neglect of the terms in dtPt[V, V], 
leading to the integrated result 



Pt[V,V]^ %,[V] 
with the decoherence rate governed by the quantity 



r = g^c2 f dt'A^'^'\t') a\t') [ dh \v{z) - V{z) ' 



(3.9) 



(3.10) 



Eq. (|3.9| ) has the form of a Gaussian functional of V{x) — V{x), whose inter- 
pretation can be estimated by evaluating it for configurations which extend over 
co-moving distances of order L and whose amplitude is Vq. For such configurations 
we have 



Pt(Vb,K)) -P^ocxp 



(3.11) 
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with a width, cr, which evolves with time according to 



= /c2 / dt' A^''-''{t')a'{t') 



(3.12) 



To proceed we assume that A evolves with time as does temperature. Then 
A(a) = Ao(ao/a) and so the width of the Gaussian density matrix acquires the time 
evolution 



where the last relationship assumes H{a) — Ho{ao/ay, as is appropriate for a 
radiation-dominated universe. 

The time-evolution of a differs qualitatively depending on whether o? > 5 or 
d < 4. If c? > 5 (corresponding to a dimcnsion-6 or higher infiaton-mattcr interaction) 
then l/cr^ saturates to g'^C2 Al'^^'^al/[{2d — 9)Ho] as a/ao grows, indicating that for 
inflaton-matter couplings this weak the Gaussian maintains an essentially fixed width 
as the universe expands. By contrast, for d < 4 we find a shrinks to zero as a/ao 
grows without bound, indicating that couplings this strong lead at late times to 
a diagonal distribution in the \V) basis, corresponding to a stochastic statistical 
ensemble of classical field configurations, V{x). 

For the two cases of most interest (a dimcnsion-5 infiaton-mattcr interaction with 
d = 4 and g = 1/M, or a dimension-4 interaction with d = 3 and g dimensionless) 
we find 



for large a/ao, where £{a) — La denotes the physical scale corresponding to the co- 
moving length L. These expressions show that a shrinks like (a/ao)~^^^ or (a/ao)'^^"^ 
in these two cases as the universe expands. 

We may also use these to estimate whether sufficient time can pass to decohere 
the modes of interest for CMB observations. To this end we evaluate the Gaussian 
distribution at the epoch of radiation-matter equality, for modes with amplitude 
Vo — Vo ~ Heq — H{aeq) and physical extent i~^{aeq) ~ i?eq, starting with an initial 
pure state at the reheat epoch (for which = Orh and Hq — H^h ~ T^^/Mp). 




(3.13) 




(3.14) 
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Using Heq ~ T^jMp, a^Jao = T^h/Teq, io = ^(aeq)(ao/aeq) = ifeqH«o/aeq) and 
Hq = ifeq(aeq/flo)^ "we have 



and 



L'H!\ fg^Al\ /Te,V 



(ifrf = 4) (3.15) 



(ifci = 3), (3.16) 



where the last approximate equahties use Aq ~ Tj-h. 

Adequate decoherence requires L^iV — VY/a'^ ^ 1, and so we see that for d = 3 
this is true at radiation-matter equahty for any V — V ^ i^eq given any reasonable 
choice for the dimensionless coupling g. It also occurs when d = A provided only that 
the scale M associated with the inflaton-matter coupling is smaller than the Planck 
scale. This shows that modes whose amplitudes differ in amplitude by as little as 
i^eq behave as if they have decohered by the time of radiation-matter equality. 

Alternatively, we can ask for what epoch does the combination L^{V — V)'^/a'^ 
become 0(1) for a mode with V — V H^q and £(aeq) ~ H~^. Using the above 
expressions we see this occurs when 

a M2 Trh 



ao M2T, 



p -^eq 
^rh 



(if d = 4) 

(ifd = 3). (3.17) 



If we demand that such modes decohere not long after the reheating epoch, such as 
when a/oo ~ for A^ not too much bigger than 1, then this requires 



and 



~ lO^Ng^ GeV (if d = 3) . (3.19) 

For reheat temperatures of this order it is not necessary to trust the solution, eq. (|3.9| ), 
for too long. 

Finally, notice that if M ~ Mp (when d = A) then we have L^{V-VY/a = 0(1) 
for y — ~ Heq regardless of the reheat temperature, and so in this case all modes 
which differ in amplitude by an amount of order ifeq would be just beginning to de- 
cohere from one another as they re-enter the horizon. It is tempting to speculate that 
this might have measurable implications for the observed temperature fluctuations 
in the CMB. 
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3.4 Decoherence During Inflation? 

We next explore the extent to which short-wavelength modes can decohere longer- 
wavelength modes during inflation itself, before a mixed environmental state has 
appeared. Our goal in this section is twofold. We flrst show why decoherence within 
this epoch cannot be reliably computed using the master-equation techniques pre- 
sented within this paper — typically because of the failure of the short-time correla- 
tion assumptions, eqs. ( |2.2[ ) and (p.4|). We then argue why it should be unlikely that 
short-wavelength modes can decohere long-wavelength ones once a reliable method 
of calculation becomes possible. 

3.4.1 Obstructions to Calculation 

Since our interest for observations is in those modes having observational implications 
for the CMB, we again ask that sector A consist of those long-wavelength modes of 
the Mukhanov fleld, v, which satisfy /c/flhc ^ H at the epoch of horizon exit. All 
other modes of the fleld v are not observed, and so can be lumped together and 
included into sector B for the purposes of CMB observations. To this end it is useful 
to divide the fleld v into a short- and long- wavelength part, v = vi + Vg, where the 
long- wavelength modes, vi, are those modes whose decoherence we wish to follow 
(sector A) and the short- wavelength modes, Vg, satisfy k/ahe S> H and provide 
part of the environment over which we trace. (We ignore those modes of v whose 
wavelengths are much longer than those which are observed, since their influence 
cannot be computed reliably using our methods of calculation.) 

Since the environment satisfles k/a^e ^ H, these modes start off in the ground 
state, |0), until a later epoch arises for which k/a ~ H, after which they start 
to become squeezed and so approach their Bunch Davies vacuum, \B). Of course, 
the environment (sector B) also includes the modes of all of the many other flelds 
which are present in the high-energy theory besides the inflaton, many of which have 
masses M > H and so for which all modes may be regarded as being 'high-energy' 
in comparison with the observable inflaton modes. For the present purposes we 
assume all of these modes to be in their Bunch-Davies vacuum state, which does not 
differ appreciably from the Minkowski vacuum to the extent that k/a ^ H and/or 
M > 

Consider, then, the following interaction lagrangian density 

Cir.t = -gB[vh,^]vf, (3.20) 

where vi denotes as above the observable inflaton modes, p is some positive-integer 
power and B is some local functional of the inflaton's short- wavelength modes, Vh, 
as well as of all other heavy flelds, ip. For instance, for the self- interactions of the 
inflaton we can use the estimate provided in Appendix ^ for the strength of its 
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lowest-dimension self-interaction, 



Ant = -^v' = + 3vf,vi + 3vHvt + vf), (3.21) 

with C, dominated (for a simple model of chaotic inflation) by the gravitational self- 
interactions, with strength 

48 aH^ 

Here H is the inflationary Hubble scale, and e is the standard inflationary slow-roll 
parameter, e = ^{V' /Vy. Since for a cubic interaction momentum conservation 
allows couplings between any three modes for which the sum (fci + k2 + k^Y — O5 
see that cubic interactions can only couple very-short to very-long wavelength modes 
if two of the modes have short wavelengths and one has long. That is, only the term 
—3 ^ v^yi which can have nontrivial matrix elements and so participate in the mixing 
of very-long and very-short wavelength modes. 

To compute the influence of sector B for the decoherence of vi we must first 
estimate the relevant correlation function, {6B{x,t)6B{x' ,t')), and see whether it 
satisfies either of the short-correlation assumptions, eqs. ( |2.2|) and (|2.4|) . We must 
do so with the environment described by a pure squeezed state, which rapidly ap- 
proaches the attractor corresponding to the Bunch-Davies vacuum, \B). For in- 
stance, starting with cubic inflaton self-interactions leads to the correlation func- 
tion {B\6v^{x,t) 6v'^{x',t')\B), while an inflaton/heavy-field interaction of the form 
-Cint = g'lpvP for a field with {B\ip\B) = would simply require the calculation of the 
Wightman function, {B\ip{x,t)ip{x',t')\B). 

The main observation to be made is that when such a correlator is computed 
using the Bunch Davies vacuum, it generically does not satisfy the condition eq. (|2.4| ). 
This is most easily seen for a heavy scalar field, ip, of mass m, for which the calculation 



has been made explicitly |Q (see also I^S])' with the result (in conformal time) 



{B\^{x,ri)tP{x',T]')\B) =[--u 



1 2^ secCvrz/) 



3 3 ^ , s 
- + z/, - - z/, 2; 1 + 



2 



2 '2 ' ' Ar^i 



4 J 167r 

(3.23) 

where s^(x, 77; a;', rj') = {r] — rj' — ieY — 6jk{x — x'Y^x — x'Y, and F[a, 6, c; z] denotes 
the standard hypergeometric function. The complex parameter, z/, is defined by 
u"^ = I — m^/if^, and becomes pure imaginary in the limit m > ^ H. This expression 
does not have the 5-correlated form of eq. (|2.4|) even in the limit m ^ H, because of 
the dependence on the separations in space and time only through the combination 
s{x,ri; x' ,7]'). The symmetries of de Sitter space ensure that this property is quite 
general for the correlations of scalar fields within their Bunch-Davies vacuum. 

One might hope to press on despite of this by returning to first principles and 
simply using an expression like eq. ( |3.23| ) in the general result, eq. ( |A.7| ), of Ap- 
pendix ^ However since this equation is derived in perturbation theory, the result 
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obtained for pA{t) in this way by integration generally cannot be trusted for times 
satisfying Vt > 0(1). For most applications this ensures that the result obtained 
cannot be used for any times of cosmological interest, reflecting the generic buildup 
of complicated correlations between sectors A and B as the sytem evolves. It is only 
because the short correlation-time assumption — Vt ^ 1, or equivalently eq. ( |2.4| ) 
— makes the time-evolution problem into a Markov process that it is possible to 
trust the solutions to eq. (^TTj) for Vt ^ 1. 

3.4.2 Why Decoherence During Inflation is Unlikely 

Although we see that it is difficult to compute decoherence effects during inflation 
within a controllable approximation, there are nevertheless strong reasons for doubt- 
ing that interactions with short-wavelength field modes during infiation can deco- 
here long-wavelength field modes. The reason hinges on the assumption that these 
short-wavelength modes are themselves in pure states, which in practice are well- 
approximated by the Bunch-Davies vacuum. This assumption is certainly a natural 
one to make for infiaton modes, given that the Bunch-Davies state is an attractor of 
the squeezing equations p . It is equally natural for the modes of other heavy fields 



(although it can be violated for some choices of scalar-field initial conditions [36|). 

The point is that for short-wavelength modes the Bunch-Davies state is very 
close to the ordinary Minkowski vacuum, and so any decoherence which is generated 
by a short- wavelength mode during infiation is also likely to be present in fiat space, 
in the absence of infiation. If it were true that having short-wavelength modes in 
their vacuum sufficed to decohere long-wavelength modes, it would be necessary to 
understand why we do not see the vacuum decohere around us all the time, and why 
quantum coherence is possible at all for the many low-energy states which have been 
studied throughout physics. Any convincing demonstration of decoherence during 
infiation by tracing out pure-state, short-wavelength modes must also explain why 
the same effect does not decohere long-wavelength modes within Minkowski space. 



4. Discussion 



In this paper we present a formalism within which it is possible to follow explicitly 
the development of decoherence within a system's long-wavelength modes due to 
their interactions with various kinds of short-wavelength physics. The formalism 
relies for its validity on there being a large hierarchy between the correlation-time 
for fiuctuations in the short-distance sector and the time-frame for evolution of the 
long-wavelength modes. 

By applying this formalism to several possible choices for the short-distance en- 
vironment we show that there is ample time during infiation for decoherence to occur 
between reheating and horizon re-entry for the modes of interest for CMB observa- 
tions. We emphasize that our calculations do not yet establish what the most im- 
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portant source of decoherence might be for inflationary perturbations. In particular, 
we cannot compute with these tools how decoherence proceeds due to environments 
whose correlation times are not much shorter than the Hubble scale, or for short- 
wavelength modes prepared within their Bunch-Davies vacuum during inflation itself, 
although, we argue that decoherence during inflation by short-wavelength modes in 
their vacuum is unlikely, given the absence of such decoherence in non-cosmological 
situations. 

Rather, our intent is merely to show that plausible sources could easily have 
sufficed to decohere the observed modes. Interestingly, if the decoherence is due to 
gravitational-strength interactions with the thermal (or other mixed) state produced 
during reheating, it could well be that the process is only just occurring for modes 
whose amplitudes differ by 0{H) as they re-enter the horizon near radiation- matter 
equality. If this should prove to be the most important source of decoherence we 
believe it would be worth exploring whether this might have observable consequences. 

We argue that the decohered system tends to a final density matrix which is 
diagonal in a basis which also diagonalizes the relevant interaction Hamiltonian. On 
grounds of locality this is typically a basis for which the long-distance field operators 
themselves are diagonal. Since such a density matrix corresponds physically to the 
establishment of a stochastic distribution of classical field configurations, it is pre- 
cisely the kind of late-time state which is implicitly assumed in standard calculations 
of the impact of inflationary fluctuations on the CMB. 



Note Added: 



While preparing this manuscript we discovered a paper, ref. ^Tj, posted to the 
arXiv by one of our collaborators on this work. It describes an analysis which overlaps 
the one presented here on the calculation of the effects of decoherence by squeezed 
states during inflation. We believe that it differs in its conclusions due to errors in 
the calculations presented there, including (but not restricted to) the application of 
the master-equation formalism outside of its domain of validity. 
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A. Interactions With a Fluctuating Environment 



In this Appendix we summarize the formahsm we use to determine the evolution 
of the density matrix for slow degrees of freedom interacting with an environment 
which contains fast fluctuations. We keep our exposition general, specializing to the 
case of interest at the end of the section. 

A.l General Formulation of the Problem 

Our goal is to provide a general master equation which dictates the evolution of 
the density matrix of a system, A, of slow degrees of freedom which interacts with 
an environment, B, of faster modes. We do not follow any of the B observables, 
and imagine that any correlations which the B variables induce into the A variables 
have correlation times which are very short compared with the time scale of the 
A-physics which is of interest. To this end we use a formalism which was developed 
for applications to physics in condensed-matter and atomic systems P^ ], 
but which has also seen service in describing neutrino evolution within complicated 



astrophysical environments 



Specifically, we imagine the system's Hilbert space to be the direct product of the 
A and B Hilbert spaces, H = Ha ^Hb, and the total Hamiltonian to be described 
by the following sum 

H = Ho + V = Ha®I + I(^Hb + V, (A.l) 

where Ha and Hb are the Hamiltonians of the subspaces A and B in the absence 
of any interactions, and V is the interaction term. In future we often suppress 
the presence of the unit operators, /. We also imagine the system starts off being 
described by a density matrix, p, for which there are no initial correlations between 
A and B: p{t = to) = ® Qb- 

In these circumstances a mean-field expansion allows practical progress, assum- 
ing that the correlation times, r, are small enough that Vt <^ 1 for the matrix 
elements of interest. The mean-field expansion must be defined in terms of the time- 
evolution operator, since only in this case do the mean-field and fluctuation parts 
not interfere with one another. That is, if we consider the time-evolution operator 
in the interaction representation, then 

p{t) = U{t)p{t = to)U*{t) (A.2) 

with 

^ = -^Vit) U{t) , (A.3) 

where V{t) = e^Hot y ^-iUot ^ 

Let us define the mean field U by 

U={U)b = Tib[pU], (A.4) 
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where the trace is only over the states in sector B. Thus, U is an operator which 
acts in Ha- Defining the fiuctuation by AU = U — U, the exact time dependence of 
any ^-sector observable, Oa, can be written as 



{OA)it) 



Tr [p{t) Oa 

TT[u{t)p{t = to)U*{t) Oa 
TTA\u{t)pA{t = to)U*{t) Oa 



(A.5) 



+Tr 



AU{t)p{t = to)AU*{t) Oy 



In the first term of the last line of this equation, pA is the reduced density matrix 

PA{t) = TTBm], (A.6) 

which suffices to follow the time-dependence of any measurement only performed in 
sector A. 

The last line of ( [A.5[ ) also shows that defining the mean in terms of U (rather 
than, say, H) guarantees that the time evolution of all observables decomposes into 
the sum of a mean evolution plus a fiuctuation about this mean, with no cross terms 
involving both U and AU. Indeed, this is the point of using U (rather than V, say) 
when making the split between mean and fiuctuating variables. 

We now compute the differential evolution of pa(^) using perturbation theory, 
with V being the perturbation. Expanding the exact evolution equation gives 



dpA 
dt 



Vit) pA{t)-pA{t) V\t) 



(A.7) 



to 



6V{t) p{t) 5V{t) + 5V{t) p{t) 5V{t) 



In these expressions 5V = V — V and V is defined in terms of U by 



so 



V{t) = {V)B-^ I dT{6V{t)6V{T))B + 0{V' 

to 



(A.8) 



(A.9) 



Notice that U need not be unitary, because probability can scatter out of sector 
A into sector B, or vice versa. Consequently, V need not be hermitian. However, an 
important special case for which V is unitary arises when the B sector involves only 
states which are too heavy to be produced by the energy available in the A sector, 
and the B sector is prepared in its vacuum state: pb = |0)bb(0|. The evolution of 
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states in sector A is guaranteed to be unitary in this case because there is no energy 
stored for release in sector B, and there is insufficient energy available in sector A to 
excite states in sector B, and to thereby complicate the density matrix in this sector. 

In the event that the system B is large enough that its state is largely unchanged 
by the interactions with A, we can evaluate the right hand side of eq. ( [A .7] ) in terms 
of pA{t) by writing p(t) ^ PAif) ® Qbi and thereby set up a differential equation 
which may be integrated to obtain the long-time behavior of PA{t)- The point of 
this formulation is that this evolution can be trusted even over very long time scales, 
provided that the correlation time of system B is sufficiently short. 

Similarly, if we write the interaction Hamiltonian in the form V{t) = Aj{t)®Bj[t) 
(with an implied sum over 'j'), and assume that the correlations in the fluctuations 
of -Bi are short-lived, 



{SBj{t)SBk{t'))B = TTB gBSB,{t)SBk{t') =W,k{mt-t' 



(A.IO) 



then we find 



V{t) = A,{t){B,{t))B - -W,k{t)A,{t)A,{t) 
and V*{t) = A,{t){Bj{t))B + W,k{t)A,{t)Ak{t) . 



(A.ll) 



These expressions use dt' f{t')5{t — t') = | f(t) and V^jkii) = VVfcj(t). The master 
evolution equation for pA{t) then is (neglecting 0(y^) terms) 



dpA 
dt 



PA, Aj {Bj)B - ^ Wjk[(^pAAjAk + AjAkPA - 2AkPAAj^ | . (A.12) 



A. 2 Microscopic Fluctuations 

We now specialize to the case where the degrees of freedom in sector B also have 
very short wavelength, in which case their infiuence can be represented in terms of 
local interactions in the spirit of effective field theories. In this case the constraints 
of locality allow considerable simplification. Suppose, then, that sector A and B 
interact through a local interaction of the form: 



V 



Aj{x,t) ® Bj{x,t) 



(A.13) 



where Aj denotes an ultra-local function of the fields describing the A degrees of 
freedom, and Bj plays a similar role for sector B. 

In terms of this interaction, the interaction potential, V, is expressible (using 
the formulae given above) in terms of correlations of the variables Bj. Using the 
assumption that these correlations in the B sector are very short compared with the 
other scales of interest, we may write 



{SBj{x, t) 6Bk{x', t'))B = Ujk{x, t) (5^(x - x') 8{t - 1') 



(A. 14) 
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where the Ujk{x,t) = U^j{x,t) are calculable local functions of position. Given this 
assumption, we have 

^ dt' {6V{t)6V{t'))B = \j d'xAjAkUjk{x,t) . (A.15) 

and so we see that V{t) = J d^xV{x,t), with 

V = A,{B,{x,t))B - '-AjAkU,k{x,t) + 0{V^) . (A.16) 

If probability and energy transfer can only flow into sector B from sector A, then 
the matrix function Ujk is non-negative definite. 

Using these expressions for the correlations of the energy density in the environ- 
ment, we can write the master equation for pAif) in the following way, accurate to 
second order in the coupling: 

—ji = i / pA,Aj {Bj)B-- / d^xUjk AjAkpA + pAAjAk -2AkpA-Aj . 

(A. 17) 

This is the result quoted as eq. ( p.5|) in the main text. 



A. 2.1 Example: Photons in a Thermal Fluid 

An explicit example of practical interest to which the key assumption, eq. (|A.14| ), 
applies is that of light propagating through an electrically neutral medium, for which 
sector A consists of the photons and sector B is the medium. Taking the photon- 
medium interaction to be given by 

£int = ev4^J^ (A. 18) 

where is the electric-current operator, shows that we may take A = eA^ and 
B = . Provided the medium is electrically neutral and carries no net currents we 
then find {J^)b = 0, and eq. ( |A.14| ) becomes 

( J^(x, t).r{x\ t'))B ~ W{x, t) 53(a; - x')5{t - t') , (A.19) 

where we follow conventional practice and use the notation 11^'^ (instead of Ujk) for 
the current correlation function. 

With these assumptions the effective photon lagrangian, eq. ( [A. 161 ) becomes 

V = -—A,A,W'{x,t) + OiV^) , (A.20) 

which shows that photon propagation through a medium is governed by the medium's 
polarization tensor, Tl^^ , with photon scattering being governed by fluctuations of 
the electric current within the medium. In the special case where the medium's 
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n 



(A.21) 



current fluctuations are due to tliermal fluctuations in a rotationally-invariant fluid 
which is in local thermal equilibrium, then they have the form 

TTi 
7^2 5'^ 

with local quantities, tti and 7r2, which are calculable in terms of the local thermo- 
dynamic variables. If the temperature is much higher than the other relevant scales 
then we expect on dimensional grounds that vr^ = kaT'^, where ki and ^2 are calcu- 
lable dimensionless constants, leading to a photon scattering rate which is of order 
r ~ e^T^. Ref. |^ shows that a similar analysis for neutrino propagation reproduces 
the standard results for neutrino scattering by a thermal bath. 

A. 3 Pure-to-Mixed Transitions 

We next identify what features of the master equation play a role in making transi- 
tions from pure states to mixed states. 

A pure state is defined by the condition that there exists a state vector, \ip), for 
which the density matrix can be written pA = This is a special case of the 

diagonal representation for pA, 



Pa 



^p„|n) (n| 



(A.22) 



for which Pn = for all n ^ ip, and = 1 for n = ip. This is equivalent to 
the basis-independent condition p\ = pA- Since pA is non-negative definite and 
initially normalized, Tr [pa] = 1, the coefficients p„ satisfy < p„ < 1 and so 
may be regarded as probabilities. In particular this implies p^ < Pn, with = p„ 
if and only if p„ = 0, 1. It follows that pA describes a pure state if and only if 
Tr [p\] = Tr [pa] = 1. 

Motivated by these observations, we next use the master evolution equation, 
( |A.12|) , to evaluate how Tr [pa] and Tr [p\] vary with time. As is simple to see, 

— Tr [pa] = i{Bj)BTT pA, Aj - - W^^Tr | (^pAAjAk + AjAkPA - 2AkPAA^ } 

= 0, (A.23) 
and the rate for the pure-to-mixed transition is 

= Tr 



dpA , dpA 



dt 



Tr 



PAAjPAAk 



dt 

Evaluating this for an initially pure state, pAito 

d_ 

dt 



PaAjA, 



■ Tr [pi] 



- 2W,fcTr 

a) (a I, then leads to 
-2Wjk{to){a\5AjSAk\a) , 



(A.24) 



(A.25) 



to 



where 5Aj 



Aj — {a\Aj\a). 
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B. Metric-Inflaton Perturbations 



Before turning to explicit choices for the environment, we first pause to describe the 
observable sector (sector A) of interest for inflationary calculations. This allows us 
to identify the modes of interest for inflationary perturbations, and to estimate the 
size of their self-couplings. 

B.l Linearized Modes 

For inflationary applications the observable system consists of those perturbations 
to the energy density which have wavelengths corresponding to those which are re- 
sponsible for the observed temperature fluctuations in the CMB spectrum. Following 
standard practice we take these to be described by the appropriate combination of 
the scalar perturbations of the metric and the modes of the inflaton field, 0, whose 
slow roll is responsible for the occurrence of inflation. Our description here draws 
extensively from the treatment in refs. |]2^, |2S[ . 

For single-field inflation the relevant dynamics is described by the action (using 
rationalized Planck units = SttG = 1) 



p 



S = — d^x 



(B.l) 



Inflation can occur if the universe should enter a phase during which the total energy 
density is dominated by a slowly-rolling homogeneous expectation value, 0(t), whose 
kinetic energy is much smaller than its potential energy: \<\P' ^ ^(0)- This can 
occur - given specific initial conditions - if the inflaton potential, ^(0), is chosen to 
be sufficiently flat. During inflation we take the space-time background to be 

ds^ = -dt^ + c?{t) dx^ dx^' , (B.2) 

where the Hubble scale, if = d/a, is approximately constant and 7jj = denotes 
the spatial metric, which for simplicity we take to be flat. 

The primordial fluctuations of relevance to the CMB within this picture are the 
quantum fluctuations of the inflaton field and of the scalar fluctuations of the metric. 
The wavelengths of these fluctuations are stretched exponentially during inflation, 
eventually to the size of the entire present-day observable universe. 

Once formed, the evolution of the fluctuation is traced by linearizing the field 
equations about the inflationary background, using = -t- 50 and g^i,v = g^j.^ + Sg^u, 
where ip and ^f^i, represent the background inflationary solution to the equations of 
motion. Under these assumptions there is a gauge for which the scalar fluctuations 
of the metric can be written in terms of a scalar field, ip, by 



{g + 6g)^, dx^dx" = a'(r/) -[1 + 2iIj{x, t])] dr/' + [1 - 2V^(x, r/)] -fijdx'dx' . (B.3) 
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We here switch to conformal time, defined by dt = a{ri) drj, and ignore vector and 
tensor metric fiuctuations since these are not yet relevant for CMB observations. 

Only one combination of the two fields 6(j) and propagates independently, and 
within the linearized approximation this combination is most efficiently identified in 
terms of the variable |3^ 



a{ri) 



(B.4) 



where primes denote differentiation with respect to r] and Ti = a' /a = Ha is the 
conformal-time Hubble scale. Expanding the action to quadratic order in v gives: 



5*2 = ^ / d'^x 



(B.5) 



where z = aif'/Ti. The utility of the variable v follows because this action has 
canonical kinetic terms, and simply describes the evolution of a free scalar field with 
a time-dependent squared mass, = —z"/z, propagating in fiat spacetime. This 
makes its evolution and quantization comparatively straightforward to understand. 

Writing the fiuctuation this way shows that the evolution of the fiuctuation 
amplitude depends crucially on the relative size of the mode wave-vector, k, and the 
'mass' term, —z"/z. In particular, this evolution changes qualitatively depending 
on whether k/a is larger or small than the Hubble scale H. What is important is 
the relative size of these two changes during the course of infiation and afterwards. 
During infiation the wavelength, Xk{t) = a{t)/k, stretches while the Hubble length, 
remains almost constant, and so the modes relevant to the later universe at 
some point leave the horizon by transiting to A > H^^ at the epoch of horizon 
exit, roughly 60 e-foldings before the end of infiation. After infiation the Hubble 
length grows faster than does A, and so modes begin to re-enter the Hubble radius 
once infiation ends (see Fig. (|1|)). After re-entry these fiuctuations are imagined 
to provide the primordial seeds which are responsible for the observed temperature 
fiuctuations in the CMB. Between horizon exit and re-entry, while k/a <C H, the 
mode amplitude satisfies v'/. ~ {z"/z)vk and so Vk{ri) oc z{ri) oc 0(77), corresponding 
to 6(f) being approximately frozen. The proportionality z oc a holds if the equation of 
state for the background geometry does not change in time, as we henceforth assume. 

The time-evolution of the quantum state of the field, v, may be understood using 
the second-quantized Hamiltonian for the action 5*2, which is 

Ha = Hfj-cc - i (^~^ j d?k{ak a-k - h.c. j , (B.6) 

where -fffree is the usual harmonic-oscillator Hamiltonian for each mode. Because the 
second term in this expression is time-dependent, it causes the energy eigenstates to 
evolve away from the adiabatic vacua, leading to the 'squeezing' for super-Hubble 
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modes which satisfy k/a <^ H . Typicalfy this evolution leads to an attractor solution, 
corresponding to the Bunch-Davies vacuum state. For sub-Hubble modes {k/a S> 
H), by contrast, the second term in (|B.6| ) is dominated by the first term, and so this 
state does not differ much from the 'Minkowski' ground state defined by afc|0) = 0. 

B.2 Coupling Strengths 

For later use we next estimate the typical strength of the self-interactions amongst 
these linearized fluctuations. It suffices for our purposes to focus on cubic self- 
interactions of the general form 

Ant = -e^', (B.7) 

for which we wish to estimate ^. Since f is a linear combination of 50 and such 
couplings can potentially arise either as infiaton or as gravitational self-interactions. 
In this section we argue that the slow-roll requirement on the infiaton potential 
ensures that gravitational interactions provide the dominant cubic contribution. 

Gravitational Self-Interactions 

An estimate the size of cubic metric self-couplings is obtained by expanding the 
Einstein-Hilbert action to cubic order in ip. Evaluating the Ricci scalar for the 



metric of eq. (|B.3|) and dropping derivatives of i}) leads to 



and so the cubic term in the Einstein-Hilbert lagrangian density, Ceh = ~\ y/—g R, 



2 

is 

24a" fn^^ 



Cs = -24 aaV = ^ H ' (^-9) 

a WJ 

where the last equality expresses the result in terms of the canonically-normalized 
variable obtained by dropping 50 from the definition oi v. ip \H/{aip'y\v = 
[H/{a^)]v. 

Specialized to the near-de Sitter geometry appropriate during inflation we have 
a" /o? = 2aH^ for constant H, related to the scalar potential by SH"^ ~ V. Fur- 
thermore, n/ip' = H/ifi ^ -3H^/V' ^ ±(2e)-^/^ where e = 1{V' /Vy is the usual 
slow-roll parameter and we use the slow-roll condition 3H(f ^ — V'. (When applied 
to y, primes here denote derivatives with respect to its argument, Eq. ( |B.9| ) then 
reduces to 



£3 = -^^a^R^^'^ = ^ 



^;^ (B.IO) 



_(2e)3/2_ 

where the upper (lower) sign applies if > (< 0). Restoring Mp = {8nG)-^/^, this 
represents a contribution to the cubic coupling which during inflation is of order 

48 aH^ 
(2e)3/2Mp 



. (B.ll) 
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Inflaton Self-Couplings 

We next compare the estimate of the cubic gravitational self-coupling just made 
with that obtained from the inflaton scalar potential. Expanding V to cubic order 
in d(j) leads to the following cubic interaction £3 = —y/^V"'{ip){6(f)Y. Since this 
depends on the third derivative, its size depends somewhat on the kind of inflationary 
potential which is used. For instance, for a simple model of chaotic inflation, defined 
by a potential of the form V = Xcj)"^, with A a dimensionless coupling constant, we 
have 

£3 = -AXa^ip{6(j)f = -4Xa^v\ (B.12) 

leading to the effective trilinear coupling constant = AXaif. 

How does this compare with the previously-computed gravitational estimate? 
Taking the ratio ^^/^g gives 



AXaip 



^g ±48a//V[(2e)3/2Mp] 



± 



12H^ \^ ^' 4 if 32 ' 

(B.13) 

where we use sigiup' = — sign(^, and simplify using the expressions SH^M^ ^ V = 
Xif^ and 2e = (V'/V)'^ = 16{Mp/ipy which are obtained using the assumed inflaton 
potential. Since the existence of a slow roll requires e <^ 1 we see that it is the 
graviton self-coupling which dominates. In what follows we therefore use a cubic 
self-interaction with strength given by the gravitational self-coupling when computing 
decoherence. 



B.3 Squeezed States 

This Appendix briefly summarizes some of the useful properties of the Bunch-Davies 
attractor, as described in the squeezed-state language. Our description here follows 
that of ref. g. 

Since the evolution of the field v whose evolution is governed by the Hamiltonian 
( p.6|) , the evolution of its ground state can be written U\0), where 



(B.14) 



U = S{n,<l>k)Riek). 
Here the operators S and R are defined by 



S{rk, = exp (e ^'^'^a^kCtk - h.c.) 



exp 



(B.15) 



and S'(rfc,$fc) defines the state's squeezing, while R{Qk) represents a rotation. 

The parameters r^, $/c and 6^ are computable functions of conformal time, 1], 
given an explicit inflationary scenario. An attractor solution for these parameters is 
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given by 0: 

1 



Tk = sinh 



-1 



2kr] 

4 \2kr] 



$fc = - tan"^ ( ) (B.16) 



Gfe = fcr? + tan-i (^^) , (B.17) 



corresponding to the Bunch-Davies vacuum of de Sitter space. 
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